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OBJECTIVE 

Determine  whether  the  ground-wave  signal  for  vlf  and  If  radio  systems  represents  the 
minimal  expected  field  strength  under  severely  disturbed  ionospheres. 

RESULTS 

Path  losses  can  markedly  exceed  those  predicted  by  normal  ground-wave  attenuation 
rates  for  propagation  under  simultaneously  occurring  conditions  of  depressed  ionospheres 
and  poorly  conducting  ground.  When  these  conditions  apply,  the  safest  way  to  evaluate  the 
vlf  and  If  system  performance  is  by  use  of  waveguide  mode  theory  rather  than  by  wave-hop 
or  ground-wave  techniques. 


RECOMMENDATION 

Because  of  the  possible  errors  and  misconceptions  which  may  arise  from  using  the 
ground-wave  computation  for  estimating  If  coverage  in  severely  disturbed  environments,  a 
complete  propagation  analysis  of  coverage  should  be  performed  by  means  of  numerical 
methods  based  on  waveguide  mode  theory. 
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1.  INTRODUCTION  AND  BACKGROUND 


Much  effort  has  been  expended  during  the  past  few  years  in  developing  survivable 
communications  for  control  of  the  strategic  forces  by  the  NCA  and  CINCS  under  nuclear 
warfare  conditions.  Both  the  Navy  and  the  Air  Force  have  developed  airborne  and  fixed 
ground-station  transmitters  which  radiate  in  the  vlf  and  Iff  10^-10^  Hz)  radio  bands  to 
provide  assured  one-way  communications  from  the  command  authorities  to  the  strategic 
forces.  The  reliance  on  long-wave  systems  arose  because  practical  transmitters  and  powers 
can  provide  near-global  (or  at  least  very-long-range)  one-way  assured  communications.  In 
addition,  long-wave  communications  were  considered  to  be  much  less  vulnerable  to  disrup¬ 
tion  by  the  effects  of  nuclear  detonations  in  the  atmosphere  than  long-range  hf  communications. 
Radio  measurements  made  during  the  1962  high-altitude  tests  in  the  Pacific  confirmed  that 
long-wave  systems  did  indeed  suffer  less  severe  propagation  effects  than  hf  systems,  but  the 
effects  were  much  more  severe  than  had  been  anticipated.  Subsequently,  researchers  were 
able  to  show  that,  if  positive  and  negative  ions  were  included  in  the  vlf  and  If  propagation 
models,  the  theoretical  predictions  were  in  line  with  observations  made  during  the  high- 
altitude  nuclear  tests. 

A  working  group  was  organized  by  the  Defense  Nuclear  Agency  (DNA)  at  the  RAND 
Corporation  in  November  1963  to  study  the  effects  of  nuclear  environments  on  If  propa¬ 
gation.  Although  details  of  the  chemistry  were  not  well  known,  it  became  apparent  to  the 
group  that  nuclear  perturbation  of  the  atmosphere  below  the  normal  D-region  could 
drastically  reduce  If  (and  vlf)  sky-wave  propagation.  The  major  conclusion  of  the  working 
group  was:  "In  any  event,  it  appears  basically  sound,  in  order  to  provide  invulnerability  to 
nuclear  burst-produced  propagation  effects,  to  design  If  systems  on  the  basis  of  the  surviv¬ 
ing  ground  wave  signal  and  normally  expected  noise  background."*  The  remarks  preceding 
the  above  quotation  are  essentially  a  distillation  of  the  introduction  and  background  from 
reference  I  and  serve  well  as  background  material  for  this  report. 

In  part  D  of  reference  I ,  Crain  illustrates  how  one  may  calculate  the  signal  amplitude 
reduction  for  a  given  path  length  and  for  several  radio  frequencies  if  the  sky  wave  were 
completely  absorbed  by  a  severe  nuclear  environment.  That  is.  the  signal  reduction  represents 
the  difference  between  signal  levels  in  a  normal  quiescent  environment  and  that  produced  by  a 
ground  wave  only.  In  presenting  this  approach.  Crain  points  out  that  it  is  approximate  and 
that  path  losses  are  apt  to  be  greater  than  given  by  the  simple  sky-wave— ground-wave  dif¬ 
ference  in  signal  level;  that  losses  increase  with  decreasing  radio  frequency  and  with  increasing 
nuclear  perturbation.  Crain  assumed  "typical"  earth  conductivity  and  typically  normal  day 
and  night  ionospheres. 

Crain's  ground-wave  attenuation  rates  were  calculated  for  free-space  conditions  above 
a  conducting  earth.  Johler“  examined  ground-wave  propagation  in  a  uniformly  ionized 
atmosphere  and  found  that  if  the  ionization  were  sufficiently  intense,  ground-wave  attenuation 


*  Crain,  C.M.,  An  Overview  Discussion  of  Propagation  Effects  of  Nuclear  Environments  on 
VLF-LF  Communications  Systems,  RAND  Corporation  Report  to  Defense  Nuclear 
Agency,  DNA  3778T.  31  August  1975. 

^  Jollier,  J.R..  Ground  Wave  Propagation  in  a  Normal  and  an  Ionized  Atmosphere,  ESSA 
Research  Laboratories  Report  ERL  1 2 1  -ITS  85.  July  1969. 
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would  be  increased.  In  a  subsequent  report,  Johler^  again  examined  the  effects  of  nuclear- 
produced  ionization  on  ground-wave  propagation  but  with  a  more  realistic  ionization  profile, 
one  which  had  a  steep  gradient  of  ionization  with  altitude.  He  again  found  sizeable  increases 
in  ground-wave  attenuation  for  large  nuclear  contamination.  Using  Johler’s  approach-^,  it 
certainly  would  be  possible  to  modify  the  path  loss  estimation  method  of  Crain  *  to  take  into 
account  ground-level  ionization  for  appropriate  nuclear  fallout  models.  However,  in  this 
report  we  examine  some  other  effects  which  increase  the  ground-wave  attenuation  of  long 
waves.  It  has  long  been  recognized  that  combinations  of  low  ground  conductivity  and  low 
ionospheric  heights  lead  to  very  large  attenuations  in  the  earth-ionosphere  waveguide. 

Fappert^  analyzed  the  effects  of  ground  conductivity  variations  on  quasi-TE  and  quasi-TM 
waveguide  modes  excited  by  both  vertical  and  horizontal  dipoles.  Calculations  for  ground 
conductivities  ranging  from  10  to  10“^  S/m  showed  that  excitation  factors  drop  rapidly  with 
decreasing  conductivity  and  that  the  dominant  TM  mode  attenuation  rates  for  10~^  S/m 
ground  conductivity  were  the  order  of  15-20  dB/Mm  as  opposed  to  a  few  dB/Mm  for  10  ^ 
S/m  ground.  The  increased  attenuation  was  much  more  severe  in  the  daytime  than  at  night. 
Field  et  al^  demonstrated  that  theoretical  computations  of  long-wave  (elf/vlf/lf)  path  losses 
over  ice  during  PC  As  were  consistent  with  experimental  measurements  of  vlf  transmissions 
over  propagation  paths  crossing  the  Greenland  ice  cap.  They  also  demonstrated  that  path 
losses  measured  in  a  scaled  physical  model  waveguide  gave  qualitative  substantiation  of  both 
the  theory  and  vlf  transmission  measurements.  Most  important,  they  conclusively  demon¬ 
strated  that  ionospheric  depression  over  low  conductivity  ground  (ice)  produces  path  losses 
far  in  excess  of  those  produced  by  the  same  ionospheric  depression  over  sea  water. 

Westerlund  et  al6  found  that  a  depressed  polar  ionosphere  from  a  proton  precipitation 
event  caused  extreme  attenuation  of  vlf  waves  passing  over  the  Greenland  ice  cap.  In  a 
subsequent  paper,  Westerlund  and  Reder ^  demonstrated  that  the  high  attenuation  of  TM 
modes  over  ice  is  a  function  of  ionospheric  height,  surface  conductivity,  and  the  dielectric 
constant  of  the  ice.  They  also  show  that  there  is  always  a  TM  mode  with  eigenangles  near 
the  Brewster  angle  for  conductivity  less  than  about  1 0  ^  S/m.  They  found  good  agreement 
between  the  theoretical  models  (using  the  NOSC  waveguide  program)  and  radio  measurements. 


Johler,  J.R.  ,  Atmospheric  Nuclear  Disturbances  with  Respect  to  the  Ground  Wave, 
ESSA  Research  Laboratories  Technical  Memorandum.  ERLTM-1TS  217,  March  1970. 

Pappert,  R.A.,  Effects  of  Elevation  and  Ground  Conductivity  on  Horizontal  Dipole 
Excitation  of  the  Earth-Ionosphere  Waveguide.  Radio  Science,  vol  5,  pp  579—590. 

March  1980. 

-*•  Field,  E.C.,  C.  Greifinger,  and  K.  Schwartz,  Transpolar  Propagation  of  Long  Radio 
Waves,  R-683-DASA,  DASA  2621.  RAND  Corporation.  March  1971. 

6-  Westerlund.  S.,  F.H.  Reder,  and  C.  Abom,  Effects  of  Polar  Cap  Absorption  Events  on 
VLF  Transmissions,  Planetary  and  Space  Science.  17.  pp  1329-1374.  19p9. 

7  Westerlund,  S.,  and  F.H.  Reder,  VLF  Radio  Signals  Propagating  over  the  Greenland 
Ice-Sheet.  Journal  of  Atmospheric  and  Terrestrial  Physics,  vol  35,  pp  1475-1491.  1973. 


In  the  following  sections,  we  look  more  closely  at  the  assumption  that  a  surviving 
ground  wave  provides  the  minimum  expected  signal  in  a  severely  disturbed  environment. 

We  briefly  discuss  the  ground-wave  model  and  indicate  weaknesses  in  the  model  for  real 
environments.  We  compare  diffraction  mode  and  waveguide  mode  attenuation  mechanisms 
as  a  function  of  ground  conductivity,  ionospheric  height,  and  frequency.  Curves  showing 
total  principal  mode  attenuation  rates  as  a  function  of  ionospheric  height,  radio  frequency, 
and  ground  conductivity  are  compared  with  ground-wave  attenuation  rates  to  indicate  the 
environmental  conditions  for  which  the  ground  wave  represents  a  reasonable,  an  unreasonably 
low,  or  an  unreasonably  high  estimate  of  the  surviving  signal  in  a  disturbed  environment. 

2.  THE  GROUND  WAVE 

The  contribution  of  ground-wave  components  to  the  total  field  is  accounted  for  by  both 
the  waveguide  mode  and  wave-hop  theories  of  EM  wave  propagation  in  the  earth-ionosphere 
waveguide.  The  mode  theory  includes  all  propagating  field  components  within  a  set  of  self- 
consistent  resonant  modes;  the  ground  wave  is  included  implicitly.  The  wave-hop  theory 
considers  the  total  field  at  some  position  in  the  guide  to  be  the  sum  of  contributions  arriving 
along  rays  connecting  transmitter  and  receiver  via  ( 1 )  reflections  off  the  earth  and  ionosphere 
and  (2)  the  ground-wave  field.  The  ground  wave  makes  an  explicit  contribution  and  is 
calculated  separately.  It  follows  from  the  wave-hop  theory  that,  if  the  ionosphere  becomes 
so  lossy  that  sky-wave  contributions  to  the  field  are  negligible,  the  ground-wave  field  would 
remain  as  the  minimum  expected  signal.  It  should  also  follow  that,  under  the  same  lossy 
conditions,  fields  calculated  by  using  waveguide  mode  theory  would  also  approach  a  minimum 
at  the  ground-wave  signal  level.  That  these  assumptions  can  be  incorrect,  particularly  under 
conditions  of  low  ground  conductivity,  is  evidenced  by  figure  1 . 

The  curve  in  figure  1  is  the  locus  of  “reflection  height”  (IT)  and  ground  conductivity 
(a)  combinations  which  produce  a  waveguide  mode  attenuation  rate  at  20  kHz  equal  to  the 
ground-wave  attenuation  rate  for  the  same  ground  conductivity.  All  combinations  of  H’  and 
a  below  the  curve  produce  waveguide  attenuation  rates  greater  than  the  ground-wave  rate; 
all  above  the  curve  produce  waveguide  attenuation  rates  less  than  the  ground-wave  attenuation 
rate.  The  attenuation  rates  considered  were  for  the  principal  waveguide  and  ground-wave  mode 
and  represent  the  total  field  attenuation  rates  well  beyond  the  horizon. 

The  waveguide  mode  parameters  here  and  throughout  this  report  are  for  ionospheric 
profiles  which  are  exponential  in  terms  of  the  conductivity  parameter  wr  where: 

wr(h)  =  ^HT=  2.5  X  105  expIMh-H’)]"1.  <>> 

as  given  by  Wait  and  Spies. ^  The  circular  plasma  frequency  Wp  is  proportional  to  the  square 
root  of  the  electron  density  and  v  is  the  electron-neutral  collision  frequency.  The  scale 
height  0  is  in  inverse  kilometers;  the  altitude  h  and  the  “reflection  height"  H'  are  in  kilometers. 
For  all  calculations,  0  =  0.25 1 5  km  * ,  taken  to  be  representative  of  the  weak  conductivity 
gradients  associated  with  the  more  lossy  disturbed  environments.  The  electron  density  profiles 
for  several  H's,  along  with  the  assumed  electron  collision  frequency  profile,  are  shown  in 
figure  2. 

®  Wait,  J.R.,  and  K.P.  Spies,  Characteristics  of  the  Earth-Ionosphere  Waveguide  for  VLF 
Radio  Waves,  US  Dept  of  Commerce.  National  Bureau  of  Standards,  Technical  Note  500. 


ELECTRON  DENSITY,  v  cm  - 

Figure  2.  Electron  collision  frequency  v  profile  and  electron  density  profiles  consistent  with 
ionospheric  conductivity  profiles  given  by  w  (h)  =  2.5  *  105  exp  [j3(h  -  H  ))s_l  where  (3  =  0.2515  km 


One  can  better  understand  the  discrepancy  between  earth-ionosphere  waveguide  and 
ground-wave  computations  of  fields  by  examining  the  ground-wave  model.  For  computation, 
the  physical  model  of  a  spherical,  finitely  conductive  earth  surrounded  by  free  space  can  be 
replaced  by  a  planar  earth  of  the  same  conductivity  under  a  continuous  refractive  atmosphere. 
The  gradient  of  refractive  index  is  such  that  a  horizontally  launched  ray  has  an  upward  curva¬ 
ture  equal  to  the  earth’s  curvature.  The  model  is  identical  to  that  often  used  for  earth- 
ionosphere  waveguides  for  the  free-space  region  below  the  ionosphere  except  that  for  the 
diffraction  model  the  free-space  region  is  infinite.  It  can  be  shown  (Budden^)  that,  for 
certain  complex  angles  of  wave  normals,  upgoing  waves  are  converted  by  an  apparent  re¬ 
flection  process  into  downgoing  waves.  At  some  angles  of  incidence,  the  downgoing  waves 
are  reflected  from  the  earth  and  become  indistinguishable  from  the  original  upgoing  waves. 
These  are,  of  course,  the  conditions  for  a  set  of  self-consistent  waveguide  modes. 

Near  the  ground,  then,  the  diffraction  modes  consist  of  both  upgoing  and  downgoing 
waves  consistent  with  the  modal  condition.  Yet  at  high  altitudes  only  upgoing  waves  are 
allowed.  This  implies  that  something  akin  to  reflection  takes  place  at  an  intermediate  height. 
The  phase  integral  method,  sometimes  used  as  a  first  approximation  for  determining  the 
modal  solutions,  involves  the  concept  of  a  “complex  reflection  height”.  ZQ.  Booker  and 
Walkinshaw*0  have  called  the  real  part  of  ZQ  the  “track  width".  The  track  width  or, 
equivalently,  “duct  height”  gives  some  physical  feel  for  the  depth  of  free  space  above  the 
ground  required  to  satisfy  the  diffraction  field  computational  model.  Figure  3  shows  the 
track  width  as  a  function  of  frequency  for  the  zero,  first,  and  second  diffraction  inodes  for  a 
standard  earth  and  perfectly  reflecting  ground. 

For  ground  of  finite  conductivity,  changes  in  earth  reflectivity  modify  the  modal 
conditions  and  consequently  the  track  width.  Figure  4  shows  the  track  width  for  the  zero- 
order  mode  as  a  function  of  radio  frequency  for  ground  conductivities  ranging  from  10  1 
S/m  to  lO--5  S/m.  The  track  widths  for  a  =  10“'*  S/m  are  nearly  identical  to  those  shown  by 
the  n  =  0  curve  of  figure  3. 

If  the  concept  of  “track  width"  has  any  real  physical  significance,  it  is  clear  that  at 
frequencies  below  about  25  kHz,  daytime  ionospheres  encroach  on  the  free-space  region 
required  by  the  ground-wave  model  for  other  than  the  zero-order  mode  (figure  3).  It  is 
equally  clear  (figure  3)  that  for  conductivities  less  than  10  ^  S/m  it  takes  but  little 
ionospheric  depression  below  normal  to  modify  the  refractivity  of  the  ground-wave  “duct" 
of  the  zero-order  mode.  The  implication  of  these  arguments  is  that,  under  conditions  lead¬ 
ing  to  complete  loss  of  sky  waves,  the  conditions  assumed  for  the  ground-wave  computation 
are  likely  not  satisfied.  Thus  the  assumption  that  the  ground-wave  field  represents  the 
minimum  possible  signal  may  be  incorrect  for  certain  combinations  of  radio  frequency, 
ground  conductivity,  and  ionospheric  height. 


Vp  Budden,  K.G..  The  Waveguide  Mode  Theory  of  Wave  Propagation.  Logos  Press.  London. 
1 96  i . 

*0-  Booker,  H.G..  and  W.  Walkinshaw,  The  Mode  Theory  of  Tropospheric  Refraction  and 
its  Relation  to  Wave-Guides  and  Diffraction,  in  Report:  Meteorological  Factors  in 
Radio  Wave  Propagation,  pp  80--I  27,  London,  Physical  Society,  1946. 
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Figure  4.  Track  width  of  the  zero-order  diffraction  mode  as  a  function  of  frequency  and  of  ground 
conductivity  a  . 


2.1  GROUND-WAVE  ATTENUATION 


Ground-wave  fields  lose  energy  by  absorption  in  the  earth  and  by  leakage  away  into 
space.  The  total  ground-wave  attenuation  rate  can  be  calculated  from  its  functional 
relationship  with  the  diffraction  mode  eigenangle.  Alternatively,  it  can  be  calculated  directly 
from  the  field  components  at  and  above  the  earth’s  surface.  Figures  5a-d  show  the  attenuations 
due  to  ground  absorption  and  to  leakage,  as  well  as  their  sum  as  a  function  of  ground  con¬ 
ductivity  for  10-kHz,  20-kHz,  40-kHz,  and  60-kHz  ground  waves,  as  calculated  from  the  field 
components  by  taking  the  height  of  the  guide  equal  to  the  track  width. 


The  curves  of  attenuation  due  to  ground  absorption  all  show  a  steady  increase  with 
decreasing  ground  conductivity  until  they  reach  a  broad  maximum  about  some  particular 
value  of  conductivity.  At  20  kHz  (figure  5b),  attenuation  by  ground  absorption  is  less  than 
1  dB/Mm  at  10~*  S/m  and  increases  to  nearly  12  dB/Mm  at  8X10“^  S/m  before  falling  off 
to  9  dB/Mm  at  10'^  S/m.  At  the  same  time,  the  attenuations  resulting  from  leakage  away 
from  the  earth  show  broad  minima  where  the  ground  absorption  losses  maximize.  In  general, 
the  leakage  losses  are  the  greater  part  of  the  total  ground-wave  attenuation  except  within  a 
frequency-dependent  band  of  surface  conductivities.  We  will  see  later  that  the  bulges  of 
maximum  ground  absorption  occur  when  the  diffraction  mode  eigenangle  is  near  the  Brewster 
angle. 


2.2  WAVEGUIDE  MODE  ATTENUATION 

Waves  propagating  in  the  earth-ionosphere  waveguide  are  attenuated  by  ohmic  “heat¬ 
ing”  of  the  ground  and  of  the  lower  ionosphere,  and  by  leakage  out  of  the  guide  along  the 
geomagnetic  field  lines.  Leakage  losses  are  negligibly  small  during  midday  or  when  the  iono¬ 
sphere  is  overionized  and  depressed.  Figures  6a-d  show  attenuation  due  to  ground  and 
ionospheric  heating  losses  for  the  principal  waveguide  mode  as  a  function  of  ground  conduc¬ 
tivity  for  several  ionospheric  heights  (H’ )  at  10,  20.  40.  and  60  kHz.  respectively.  The  total 
ground-wave  attenuation  rate  is  also  included  in  each  of  the  figures. 

There  are  some  striking  differences  in  the  attenuation  due  to  ground  heating  between 
the  ground-wave  and  waveguide  modes.  For  both  models,  the  attenuation  from  ground 
absorption  has  a  broad  maximum  centered  about  those  conductivities  where  the  waveguide 
and  diffraction  mode  eigenangles  approach  the  surface  complex  Brewster  angle.  However, 
the  waveguide  modal  eigenangles  approach  the  Brewster  angles  at  somewhat  lower  conduc¬ 
tivities  than  do  the  diffraction  mode  eigenangles,  particularly  at  the  lower  frequencies.  The 
most  important  difference  obviously  results  from  the  presence  of  an  upper  boundary  of  the 
waveguide.  Lowering  the  ionosphere  not  only  increases  ionospheric  absorption,  but  also 
greatly  increases  ground  absorption.  At  the  lower  conductivities,  even  if  the  assumption 
were  true  that  sky-wave  absorption  is  identically  equal  to  the  ground-wave  leakage  losses,  the 
enhanced  ground  losses  under  the  depressed-ionosphere  can  materially  exceed  ground-wave 
surface  absorption.  In  the  next  section  we  will,  within  the  constraints  of  the  model  iono¬ 
sphere  show  the  conditions  undeT  which  waveguide  attenuation  rates  will  exceed  ground- 
wave  attenuation  rates. 


12 


tmt* 


Figure  6  a-d.  Attenuation  rate  of  the  principal  waveguide  mode  resulting  from  ground  absorption  (solid 
curve)  and  resulting  from  ionospheric  absorption  (dashed  curve)  as  a  function  of  ground  conductivity 
for  a  number  of  reflection  heights  H'.  The  broken  curve  is  the  ground-wave  attenuation  as  a  function  of 
surface  conductivity. 
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3.  DISCUSSION  OF  RESULTS 


3.1  COMPARISON  OF  WAVEGUIDE  AND  GROUND  WAVE 


Waveguide  mode  attenuation  rates  were  calculated  over  the  frequency  band  of  10 
to  60  GHz  by  using  the  NOSC -developed  “MODESRCH”^  *  program.  Ground-wave  attenu¬ 
ation  rates  were  calculated  by  using  an  adaptation  of  a  computer  program  developed  by  Berry 
and  Herman.  *  “  Both  programs  also  were  used  to  calculate  attenuation  rates  for  a  number 
of  ground  conductivities,  but  the  waveguide  calculations  were  also  parametric  in  terms  of  the 
“reflection  height”  for  exponential  ionospheric  profiles.  The  results  of  the  calculations  are 
shown  in  figures  7a  through  7g,  where  the  curves  in  each  figure  are  for  a  specific  ground  con¬ 
ductivity.  The  conductivities  range  from  10'*  S/m  (figure  7a)  to  10~^  S/m  (figure  7g). 

The  attenuation  rate  for  a  mode  which  is  characterized  by  a  complex  eigenangle  at  the 
ground  equal  to  the  Brewster  angle  is  given  by  the  curve  labeled  6q. 

Ionospheric  anisotropy  was  included  in  the  waveguide  mode  calculations.  All  results 
are  for  a  dip  angle  of  51°,  a  magnetic  flux  density  of  4.31  X  10“^  webers  per  square  meter, 
and  a  geomagnetic  azimuth  of  propagation  of  58.5°.  However,  because  most  of  the  calcu¬ 
lations  are  for  depressed  daytime  ionospheres,  the  magnetic  effects  are  negligible  and  the 
results  are  representative  of  isotropic  ionospheres. 

From  figures  7a-g,  it  is  clear  that  path  losses  under  depressed  ionospheres  can  materi¬ 
ally  exceed  the  theoretical  ground-wave  losses,  particularly  at  the  low  end  of  the  vlf-lf  spectrum. 
However,  for  surface  conductivities  equal  to  or  greater  than  about  10  ^  S/m,  the  ionosphere 
must  be  extraordinarily  depressed  before  waveguide  attenuation  rates  exceed  ground-wave 
rates,  particularly  at  If. 

At  conductivities  below  10~-*  S/m,  the  waveguide  attenuation  rates  begin  to  exceed 
ground-wave  rates  for  moderately  depressed  ionospheres.  In  the  3  X  10  ^  to  10~^  S/m 
conductivity  range,  the  phenomenon  pointed  out  by  Westerlund  and  Reder^  becomes 
evident.  That  is,  the  maximum  attenuation  across  the  frequency  band  is  closely  associated 
with  the  modal  eigenangles  near  the  Brewster  angle. 


The  data  described  by  figures  7a-g  are  summarized  from  a  different  perspective  in 
figure  8,  where  those  combinations  of  reflection  height  H'  and  surface  conductivity  o.  which 
lie  above  each  frequency  curve,  produce  a  smaller  attenuation  rate  for  the  principal  TM  mode 
than  the  ground  wave.  For  combinations  of  H'  and  a  lying  below  each  frequency  curve,  the 
principal  earth-ionosphere  waveguide  mode  attenuation  rate  is  greater  than  that  of  the  ground 
wave. 


Morfitt,  D.G.,  and  C.H.  Shellman,  ‘MODESRCH’,  an  Improved  Computer  Program 
for  Obtaining  ELF/VLF/LF  Mode  Constants  in  an  Earth-Ionosphere  Waveguide. 
Naval  Electronics  Laboratory  Center  Interim  Report  No.  77T  for  DNA.  1  October 
1976. 

Berry,  L.A.,  and  J.E.  Herman,  A  Wave  Hop  Propagation  Program  for  an  Anisotropic 
Ionosphere,  Telecommunications  Research  Report  OT/ITS  RR1  I,  US  Department 
of  Commerce,  Office  of  Telecommunications,  Institute  of  Telecommunications 
Sciences,  Boulder,  CO,  April  1971. 
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Figure  7  a-g.  Attenuation  rate  of  the  principal  waveguide  mode  as  a  function  of  frequency  for  a  numhci 
of  reflection  heights  H’  and  for  the  ground  conductivity  o  specified  on  each  figure.  The  ground-wave 
attenuation  rate  (curve  GW)  and  the  attenuation  rate  of  a  mode  characterized  by  eigenangles  equal  to 
the  complex  Brewster  angle  (curve  0g)  are  alto  shown  for  each  ground  conductivity 


a.  dB/Mm  o,  dB/Mm 


a.  dB/Mm 


FREQUENCY,  kHz 

FIGURE  71' 


Of,  dB/Mm 


ire  8.  Combinations  of  reflection  height  H'  and  surface  conductivity  a  which  produce  a  principal  TM  waveguide  mode  attenuation  rate  equal 
he  ground-wave  attenuation  rate.  For  those  combinations  of  H'  and  a  lying  above  each  of  the  frequency  curves,  the  TM  waveguide  mode 
nuation  is  less  than  the  ground-wave  attenuation  rate  for  that  frequency  and  ground  conductivity.  For  combinations  below  the  curve,  the 
mode  attenuation  rate  is  greater  than  the  ground-wave  attenuation  rate. 


For  the  particular  ionospheric  conductivity  gradient  (3  used  tor  this  study  and  lor 
ground  conductivity  a  less  than  about  10  ^  S/m.  ionospheric  depressions  typical  of  polar 
cap  events  or  nuclear  disturbances  can  produce  path  losses  greater  than  ground-wave  path 
losses.  Although  the  conductivity  gradient  (3  =  0.25 15  km  *  was  selected  after  an  examin¬ 
ation  of  a  number  of  conductivity  profiles  for  variously  disturbed  ionospheres,  it  should  be 
pointed  out  that  the  model  becomes  less  appropriate  for  the  more  normal  reflection  heights. 
For  example,  figure  7g  shows  waveguide  mode  attenuation  rates  at  10  kHz  for  reflect  ion 
heights  up  to  95  km.  From  the  data,  one  deduces  that  for  10  **  S/m  ground  conductivity 
and  reflection  heights  below  90  km.  waveguide  path  losses  exceed  ground-wave  path  losses. 
But  it  is  hard  to  imagine  the  circumstances  which  would  produce  a  scale  height  (3  =  0.25  I  5 
km  '  for  altitudes  much  greater  than  70  km.  Thus  the  results  shown  in  figures  7  and  N  are 
most  useful  for  locating  possible  problem  areas  but  not  their  precise  boundaries. 

4.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  assumption  that  If  communications  systems  designed  to  operate  at  ground-wave 
signal  levels  will  maintain  connectivity  in  a  severe  nuclear  environment  is  safe  only  for  pro¬ 
pagation  over  relatively  highly  conductive  ground.  At  the  lower  vlf  frequencies,  the  assump¬ 
tion  is  never  safe.  Lf  systems  designed  to  operate  on  ground-wave  propagation  over  highly 
conductive  ground  should  have  a  considerable  system  margin.  For  If  propagation  over  poorly 
conducting  ground  or  fresh  water  ice.  losses  may  considerably  exceed  ground-wave  losses 
when  the  ionosphere  is  moderately  depressed. 

When  the  DNA  working  group  recommended  that  If  systems  be  designed  on  the 
basis  of  a  surviving  ground  wave,  it  was  recognized  that  methods  for  calculating  If  propagation 
and  the  propagation  environment  were  inadequate.  Today,  although  we  still  have  questions 
regarding  some  of  the  air  chemistry,  there  is  little  doubt  that  severe  disturbances  can  be 
modeled  with  good  confidence.  In  addition,  we  no  longer  rely  on  “wave-hop''  methods 
(with  the  attendant  inaccuracies  for  severely  disturbed  ionospheres)  for  calculating  If  pro¬ 
pagation  parameters.  Waveguide  mode  propagation  prediction  methods  work  well  to  at  least 
100  kHz.  Because  of  the  possible  errors  and  misconceptions  which  may  arise  from  using  the 
ground-wave  computation  for  estimating  If  coverage  in  severely  disturbed  environments,  it 
seems  far  safer  to  do  a  complete  propagation  analysis  of  coverage  by  means  of  numerical 
methods  based  on  waveguide  mode  theory. 
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CODE  402 
CODE  420 
CODE  421 
CODE  461 
CODE  464 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CENTER 
4301  SUITLAND  RD  BLDG  5 
WASHINGTON.  DC  20390 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 
CODE  81  (HD  SMITH) 

CODE  532  (31 

CODE  532  (WILLIAM  F  MOLER) 

COMMANDING  OFFICER 
NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  DC  20375 

CODE  5410  (JOHN  DAVIS) 

CODE  7701  (JACK  0  BROWN) 

CODE  5461  TRANS  IONO  PROP 
COOE  5465  PROP  APPLICATIONS 
CODE  5460  ELECTROMAG  PROP  BR 
CODE  2600  TECH  LIBRARY  12) 

OFFICER  IN  CHARGE 
WHITE  OAK  LABORATORY 
NAVAL  SURFACE  WEAPONS  CENTER 
SILVER  SPRING,  MO  20910 

CODE  WA501  NAVY  NUC  PRGMS  OFF 
CODE  WX21  TECH  LIBRARY 


COMMANDER 

NAVAL  TE LECOMMUNCATIONS  COMMAND 
NAVTELCOM  HEADQUARTERS 
4401  MASSACHUSETTS  AVE,  NW 
WASHINGTON,  DC  20390 
CODE  24C 

COMMANDING  OFFICER 
NAVY  UNDERWATER  SOUND  LABORATORY 
FORT  TRUMBULL 
NEW  LONDON.  CT  06320 
PETER  BANNISTER 
DA  MILLER 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
NAVY  DEPARTMENT 
WASHINGTON,  DC  20376 
NSP-2141 

DEPARTMENT  OF  THE  AIR  FORCE 
CBmmanBER  ' 

ADC/DC 

ENT  AFB,  CO  80912 
DC  (MR  LONG) 

COMMANDER 
ADCOM/XPD 
ENT  AFB,  CO  80912 
XPQDQ 
XP 

AF  GEOPHYSICS  LABORATORY.  AFSC 
HASCOM  AFB.  MA  01731 
CRU  <S  HOROWITZ) 

PHP  (JULES  AARONS) 

OPR  (JAMES  C  ULWICK) 

OPR  (ALVA  T  STAIR) 

SUOL  (RESEARCH  LIBRARY)  (2) 

AF  WEAPONS  LABORATORY,  AFSC 
KIRTLAND  AFB,  NM  87117 
SUL  (2) 

SAS  (JOHN  M  KAMM) 

DYC  (CAPT  L  WITTWER) 

AFTAC 

PATRICK  AFB,  FL  32925 
TN 
TD-3 
TD-5 

TF/MAJ  WILEY 

AIR  FORCE  AVIONICS  LABORATORY,  AFSC 
WRIGHT  PATTERSON  AFB,  OH  45433 
A  AD 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION.  AFSC 
WRIGHT  PATTERSON  AFB,  OH  45433 
ETD  BL  BALLARD 

HQ  USAF/RD 
WASHINGTON,  DC  20330 
RDQ 


HEADQUARTERS 

NORTH  AMERICAN  AIR  DEFENSE  COMMAND 
1500  EAST  BOULDER 
COLORADO  SPRINGS,  CO  80912 
CHIEF  SCIENTIST 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFrISS  AFB,  NY  13440 
EMTLD  DOC  LIBRARY 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 
EEP  JOHN  RASMUSSEN 

SAMSO/MN 

NORTON  AFB,  CA  92409 

MINUTEMAN  (NMML  LTC  KENNEDY) 


COMMANDER  IN  CHIEF 
STRATEGIC  AIR  COMMAND 
OFFUTT  AFB,  NB68113 
NRT 

XPFS  (MAJ  BRIAN  G  STEPHANI 
DOK  (  CHIEF  SCIENTIST) 

US  ENERGY  RESEARCH  AND  DEV  ADMIN 
DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
PO  BOX  5400 

ALBUQUERQUE,  NM  87115 

DOC  CON  FOR  D  W  SHERWOOD 

DIVISION  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON.  DC  20545 

DOC  CON  FOR  DONALD  I  GALE 

LAWRENCE  LIVERMORE  LABORATORY 
PO  BOX  808 

LIVERMORE.  CA  94550 

GLENN  C  WERTH  L-216 
TECH  INFO  DEPT  L-3 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 

PO  BOX  1663 

LOS  ALAMOS.  NM  87545 

DOC  CON  FOR  T  F  TASCHEK 
DOC  CON  FOR  D  R  WESTERVELT 
DOC  CON  FOR  P  W  KEATON 
DOC  CON  FOR  J  H  COON 

SANOIA  LABORATORIES 
LIVERMORE  LABORATORY 
PO  BOX  969 

LIVERMORE,  CA  94550 

OOC  CON  FOR  B  E  MURPHEY 
DOC  CON  FOR  T  B  COOK  ORG  8000 


SANDIA  LABORATORIES 
PO  BOX  5800 

ALBUQUERQUE.  NM  87115 

DOC  CON  FOR  SPACE  PROJ  DIV 

DOC  CON  FOR  A  D  THORNBROUGH,  ORG  1245 

OOC  CON  FOR  W  C  MYRA 

DOC  CON  FOR  3141  SANDIA  RPT  COLL 

OTHER  GOVERNMENT 
DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  DC  20234 
RAYMOND  T  MOORE 

DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TELECOMMUNICATIONS 
INSTITUTE  FOR  TELCOM  SCIENCE 
BOULDER.  CO  80302 
WILLIAM  F  UTLAUT 
L  A  BERRY 
A  GLENN  JEAN 
0  D  CROMBIE 
J  R  WAIT 

DEPARTMENT  OF  TRANSPORTATION 
OFFICE  OF  THE  SECRETARY 
TAD-44.1,  ROOM  10402  B 
400  7TH  STREET,  SW 
WASHINGTON,  DC  20590 
R  L  LEWIS 
R  H  DOHERTY 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 
AEROSPACE  CORPORATION 
PO  BOX  92957 
LOS  ANGELES.  CA  90009 
IRVING  M  GARFUNKEL 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  RD 
BURLINGTON,  MA  01803 
RADIO  SCIENCES 

THE  BOEING  COMPANY 
PO  BOX  3707 
SEATTLE,  WA  98124 
GLENN  A  HALL 
J  F  KENNEY 


UNIVERSITY  OF  CALIFORNIA 
AT  SAN  DIEGO 

MARINE  PHYSICAL  LAB  OF  THE 
SCRIPPS  INSTITUTE  OF  OCEANOGRAPHY 
SAN  DIEGO,  CA  92132 
HENRY  G  BOOKER 

COMPUTER  SCIENCES  CORPORATION 
PO  BOX  530 

6565  ARLINGTON  BLVD 
FALLS  CHUCH.VA  22046 
D  BLUMBERG 

UNIVERSITY  OF  DENVER 
COLORADO  SEMINARY 
OENVER  RESEARCH  INSTITUTE 
PO  BOX  10127 
DENVER, CO  80210 
DONALD  DUBBERT 
HERBERT  REND 
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ESL,  INC 
496  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
PO  BOX  8555 

PHILADELPHIA,  PA  19101 

SPACE  SCIENCE  LAB  (MH  BORTNER) 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET 
PO  DRAWER  QQ 
SANTA  BARBARA  CA  93102 
B  GAMBILL 

DASIAC  (21 

■“  DON  CHANDLER 
WARREN  S  KNAPP 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
T  N  DAVIS 
NEAL  BROWN 
TECHNICAL  LABORATORY 

GTE  SYLVANIA,  INC 
ELECTRONICS  SYSTEMS  GRP 
EASTERN  DIVISION 
77  A  STREET 
NEEDHAM,  MA  02194 
MARSHAL  CROSS 

I  IT  research  institute 
10  WEST  35TH  STREET 
CHICAGO,  IL  60616 

TECHNICAL  LIBRARY 

UNIVERSITY  OF  ILLINOIS 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

URBANA,  IL  61803 

AERONOMY  LABORATORY  (2) 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL.  MD  20810 
J  NEWLAND 
PT  KOMISKE 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
3261  HANOVER  STREET 
PALO  ALTO,  CA  94304 
E  E  GAINES 
W  L  IMHOF  D/52-12 
J  B  REAGAN  D652  12 
R  G  JOHNSON  D/52  12 

LOWEL  RESEARCH  FOUNDATION 
450  AIKEN  STREET 
LOWELL.  MA  01864 
DR  BIBL 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 
PO  BOX  73 

LEXINGTON,  MA  02173 
DAVE  WHITE 
J  H  PANNE LL  L  246 
D  M  TOWLE 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
R  HENDRICK 
F  FAJEN 
M  SCHEIBE 
J  GILBERT 
C  L  LONGMIRE 

MITRE  CORPORATION 
PO  BOX  208 
BEDFORD.  MA  01730 
G  HARDING 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD 
SANTA  MONICA,  CA  90404 
E  C  FIELD,  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERIC  RESEARCH  LABORATORY 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

IONOSPHERIC  RSCH  LAB  121 

R&D  ASSOCIATES 
PO  BOX  9695 

MARINA  DEL  REY,  CA  90291 
FORREST  GILMORE 
WILLIAM  J  KARZAS 
PHYLLIS  GREIFINGER 
CARL  GREIFINGER 
A  ORY 

BRYAN  GABBARD 
R  P  TURCO 
SAUL  ALTSCHULER 

RAND  CORPORATION 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

TECHNICAL  LIBRARY  121 

CULLEN  CRAIN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
DONALD  NEILSON 
GEORGE  CARPENTER 
W  G  CHETNUT 
J  R  PETERSON 
GARY  PRICE 

STANFORD  UNIVERSITY 
RADIO  SCIENCE  LABORATORY 
STANFORD,  CA  94305 
R  A  HELLIWELL 
FRASER  SMITH 
J  KATSUFRAKIS 

TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
DIANA  DEE 

CALIFORNIA  INSTITUTE  OF  TECHNOl  OGY 
JET  PROPULSION  LABORATORY 
4800  OAK  GROVE  DRIVE 
PASADENA,  CA  91103 
ERNEST  K  SMITH 
(MAIL  CODE  144  B13I 
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